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A complete, unitied description is given ot the design, 
implementation and use of a family of very fast and efficient 
large scale minimum-cost (primal simplex) network programs. 

The class ot capacitated generalized transshipment problems 
solved includes the capacitated and uncapacitated generalized 
transportation problems and the continuous generalized assign- 
ment problem, as well as the pure network flow models which are 
specializations ot these problems. These formulations are used 
tor a large number of diverse applications to determine how (or 
at what rate) flows through the arcs of a network can minimize 
total shipment costs. A generalized network problem can also 
be viewed as a linear program with at most two non-zero entries 
in each column ot the constraint matrix; this property is ex- 
ploited in the mathematical presentation with special emphasis 
on data structures for basis representation, basis manipulation, 
and pricing mechanisms. a literature review accompanies compu- 
tational testing ot promising ideas, and extensive experimenta- 
tion is reported which has produced GENNBT, an extreioely effi- 
cient family of generalized network systems. 



1 . iNTKOUUCTiUN 

This paper reports the development ot a larye-scaie primal 
network code tor solving capacitated generalized transshipment 
problems. The capacitated generalized transshipment problem is 
the most general of the minimum cost flow models in continuous 
variables, which include the capacitated and uncapacitated trans- 
portation problems and the continuous generalized assignment problem 
as well as the pure network specializations ot these problems. 

These models are used tor a large number ot diverse ap^jlicat ions 
that include transportation ot goods, design ot reservoir, com- 
munications, and pipeline systems, assignment ot personnel ana 
machinery to jobs, bid evaluation; currency exchange and cash 
management, production, sales and inventory planning; and many 
others. For further discussion ot these applications see survey 
articles such as Bradley [bj, or Glover et ai., [17, 18], and 

textbooks (as well as their cited references) such as Uantzig 
[13], Jensen and Barnes [23], and Kennington and Helgason [27]. 

The capacitated generalized transshipment model and its 
specializations are minimum-cost network flow problems. The 
goal is to determine how (or at what rate) flows through the 
arcs ot a network can minimize shipment costs. The network is a 
directed graph, G , defined by a set ot nodes, N , and a set ot 
arcs, A , with ordered pairs ot nodes (tail, head) as elements 
indexed by k: (tj^,t|_,). For each arc there is a shipping cost 

per unit flow, c^^ , a minimum allowable flow (or lower bound), 

, and a maximum allowable tlow (or upper bound, or capacity), 
u 

k 



In addition, there are coefficients (or multipliers, or gains, 
or losses), a. and b. which can change the magnitude of (or 

K “"K 

amplify, or attenuate) each unit of flow respectively entering 
and leaving arc k . Each node is either a supply node where 
units of the good enter the network, a demand node where units 
leave, or a transshipment node. The problem is to minimize 
total costs with flows, , that satisfy the associated lower 

and upper bounds and preserve the conservation of flow at each 
node : 



(GN) 



MIN I 



keA 



c, X, 
k k 



s . t . 



I 

keA 

with t|^ = i 



^k^k 



^ -k^k 
keA 

with t, =i 
k 



X, = 



bi , 



i e N 



< X, < u, , keA, 

— k — k — 

is a supply node; 
i is a demand node; 

Note that any linear program with at most two nonzero coefficients 
associated with each variable is a generalized network (GN). for- 
mulation (Gn) is a generalized transshipment model (GT) if all 
are +1, in which case the corresponding coefficient is called the 
multiplier m. = -b. . For purposes of exposition, we will address 



where : 



bi = 



supply if i 
- demand it 
U otherwise. 
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(GT) since assuming the existence of a finite upper bound on each 
variable it is possible to transform the (GN) coet t icients--by 
scaling or by reflecting variables with respect to their upper 
bounds--so tliat one coefficient is +i tor each variable. The other 
coefficient for each variable then becomes the multiplier value, 
m. , and the generalized transshipment (GT) network flow interpre- 
tation results with a node tor each constraint and a directed arc 
tor each variable. It a variable has two nonzero coefficients, its 
arc is directed away from the node corresponding to the constraint 
with the +1 coefficient; a variable with just one nonzero coefficient 
(±1) corresponds to an arc forming a self-cycle, leaving and 
returning to the same. node. 

Some generalized networks (GN), such as those obtained by 
relaxing integer restrictions on flow variables, cannot be scaled 
conveniently to (GT). These (GN) are accommodated by obvious minor 
modifications in the following (GT) presentation. We have developed 
codes tor solving (GN) and specialized them tor (GT) problems. 

Generalized networks can be solved as linear programming 
problems, but contemporary commercial linear programming sys- 
tems consume much more computer time and data storage region 
than special purpose network codes. Indeed, the advantage of 
network codes is so pronounced that it is even worthwhile to 
develop special linear programming procedures to exploit intrin- 
sic network structure found embedded within more general models 
(e.g., Kennington and Helgason [27], McBride [3U], Brown and Graves 
[8], and Brown and McBride [9]). Brown and Wright [11] and brown, 
McBride and Wood [10] show that many real-life linear programs 
contain a large embedded generalized network structure. 
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These models are widely used because they accurately de- 
scribe a large variety of important applications. Generalized 
networks not only directly represent gains with > 1 (e.g. 

interest return on investment, heat gain, etc.) and losses with 
0 < m, <1 (e.g., evaporation, voltage drop, attrition, etc.) 

-rk 

in the flows, but also admit conversions of units for these flows 
(e.g., maqhine time to output pieces, lira to yen, etc.). In 
addition, m. < 0 represents situations without obvious physical 
flow interpretation (e.g., flows which "enter the head" of arc k 
in proportion of those entering the tail), but which nonethe- 

less provide valuable modelling tools. There has been continuing 
growth of interest in network models because efficient computer 
programs have made possible the reliable, economic solution of 
problems with more variables than virtually any other optimiza- 
tion technique (e.g., the pure network system, GNET [6], has been 
installed at hundreds of sites worldwide and is now cited as a 
routine research tool). Perhaps most important, networks are 
readily accepted by nonanalysts and are consequently extremely 
popular operations research models. 

Although several papers have been written in this general 
area, and significant computational breakthroughs have been 
reported, there has not previously been a single, unified de- 
scription of a complete implementation, nor have "new generation" 
computer programs been made generally available to the academic 
community. Here we report the research and computational experi- 
ments which have produced GENiMET, an extremely efficient family 
of network optimization systems. GENNET exploits pure network 
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structure embeddea in generalized networks, and specializes in- 
trinsically to GNET [b], when both systems are applied to pure 
networks, (the floating point arithmetic in) GENNeT requires 
about lb percent more time than GNET. An important objective of 
this paper is to make these new approaches easily accessible to 
a wide audience via a clear exposition and concrete examples of 
efficient FORTRAN programs. Further, the availability of the 
(GT) and (GN) computer programs will now make it possible tor 
other investigators to reproduce and extend our experimental 
results . 

Bradley, Brown and Graves IbJ trace the historical develop- 
ments leading to contemporary primal simplex pure network algo- 
rithms, their supporting data structures, and efficient imple- 
mentations such as GNET. For other sub-classes of generalized 
networks, algorithms have been reported by Jewell [2b] , Eisemann 
[14], Maurras 129J, Glover, Hultz, Klingman and Btutz [17,la,ly], 
Balachandran [2], and Jensen and Bhaumik I24j. An efficient 
algorithm for large generalized network problems has been de- 
veloped by Glover, Klingman, Hultz, Stutz, Karney and Elam 
( lb , 1 7 , 1 B , 21 , 2 2 ] . However, their contributions are scattered 
among the papers referenced; Kennington and Helgason [27J and 
Jensen and Barnes [23J provide textbook descriptions of compu- 
tations apparently gleaned from these papers, providing an ade- 
quate treatment of the graphical algorithm with more computational 
advice than the seminal presentation by bantzig [IJ] but tew 
details of efficient basis updating. 



b 



the; approach 



2 . 

Our approach continues with generalized networks in pre- 
cisely the philosophical vein ot the pure network exposition 
of Bradley, Brown and Graves: we seek data structures and algo- 

rithms that yield efficient implementations without abandoning 
the flexibility of a general large-scale mathematical program- 
ming perspective [6,p. 3 ff]. We introduce few of the details of 
the general bounded-variable simplex algorithm, and we repeat 
little of the underlying pure network material; the assiduous 
reader might well review the prior paper for which this is an 
intimate companion. 

We continue with a brief description of the algebraic 
specialization of the simplex method for generalized networks. 
Specific design decisions and experiments carried out with GENNBT 
are described, including computational tests ot alternate ap- 
proaches. Some extensions of GENNET are presented to further 
exploit special problem structure. 
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3. 



GENERALIZED NETWORK SPECIALIZATION 



Efficient primal simplex specialization to the generalized 
network case depends upon the well-known result that any 
generalized network basis can be put in ne^^rly (upper) tri- 
angular form by simple permutation of rows and columns. This 
inherent near triangularity can be exploited by direct solution 
of the simplex equations with mod i f ied forward, or back substi- 
tution. Fortuitously, this basis structure also leads to ex- 
tremely fast solution updates orchestrated in concert with 
efficient dynamic reorganization of each new basis. 

Theorem (e.q., Dantzig [13]) Any basis b extracted from a 
generalized network problem can be put in the form (1) by 
rearranging rows and columns . 




B 



I 



( 1 ) 






I 

where each square submatrix component B is either upper 
triangular or nearly upper triangular with only one element 
below the diagonal . 
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Proof. Our proof is constructive. We know that the basis has at 



least one nonzero entry in each row, and one, or two, 
entries in each column. Define a pairing as the associ- 
ation of a row with a column sharing an entry in B , a 
deferral as a temporary deletion of a pair from consider- 
ation, and an ass iy nment as the fixing of a pair on the 
diagonal in the ordering of the rearranged sequence, 
followed by the assignment of all deferred pairs which 
have a column with an entry in an assigned row. 

Step 1) Defer singleton rows . Locate a row with one entry, pair 
the row with the column of the entry and defer the pair. 
Repeat Step .1 until no row remains with one entry. 

Comment: Step 1 reduces B to a submatrix with exactly two 

entries in each row and in each column . To see this, 
note that each column can have at most two entries, or 
2m entries for m columns. Bach row has at least two 
entries. Suppose that some row has more than two 
entries; then at least 2m + 1 entries exist, leading to 
a contradiction. Thus, exactly 2m entries remain; 
consequently, each column has exactly two entries, as 
does each row. 

Comment: Step 1 defers an upper triangular set of pairs. To 

see this, diagonalize the pairs in reverse of the order 
of their deferral, and place this sequence at the end 
of the rearrangement of B . 

Comment: A sequence of assignments in the rearrangement begins. 
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step 2) Assign a Component . 



2.1) It a row remains which is neither assi’^neo nor deterred, 

begin a near triangular component: pa i r the row and a 

column and assign this pair next in the rearranged 
sequence. otherwise, go to Step 2.3. 

2.2) Apply Step 1 and assign each newly deterred pair next in 
the sequence. Go to Step 2.4. 

2.3) Begin a triangular component: ass ign the most recently 

deterred pair next in the sequence. 

2.4) Repeat Step 2 until all rows and columns are assigned. 

Comment: It Step 2.1 assigns a pair to a component, then the 
component is nearly triangular with one element below 
the diagonal in the first column. otherwise, the 
component is triangular. 

Consider the generalized network problem given in figure 1, 
a generalized transshipment problem with 2 sources, lu sinks, 

15 nodes, and 3U arcs. We will use this problem to illustrate 
concepts and efficient solution methods. In this problem the 
multipliers are all positive. for arc k directed trom node 
i to node j , if flow X|^ leaves node i , then 
arrives at node j . V\/hen > 1 the amount arriving at node 

j is greater than the amount leaving node i . This would be 
the case in cash flow problems when the arc corresponds to an 
investment ana mj^ = (1 + ) with the rate of return. 

When mj^ < 1 then the amount arriving at node j is less than 
the amount leaving node i . here the loss could correspond to 



y 



evaporation, taxation, transmission loss, brokeraye tees, seepage 
or deterioration. Pure network arcs are indicated by 1 

(and may be even more abundant in real problems than in our 
example,). For clarity, minimum allowable flows are zero, and 



maximum 


allowable 


flows are 


not 


specif ied . 
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Figure 1. A Single Commodity Generalized Transshipment 

Problem (GT) 
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b'igure 2 shows a basis for the problem introduced in 
Figure 1. (In this simple case, there is only one component: 

p = 1 . ) 

A unique subgraph partition ot b denoted corre- 
sponds to ts . Let = {a^|a^ is an arc associated with b-* , 

a column ot b} , then = [N,Aj^J denotes the directed graph 

0 

associated with the basis B . To each submatrix b ot b 

£ £ £ 

there corresponds a component ot G^^ denoted by G = [N ,a ] . 

£ £ £ 

N is the set ot nodes corresponding to the rows ot b and A 

£ 

is the set ot arcs corresponding to the columns ot B . It is 



known 


(e .g . , 


[13] : 


) that 


IS 


either a rooted 


tree or 


a one- 


tree 


(a tree 


with 


an additional 


arc 


forming one 


cycle ) . 


It G^ 


is a 


rooted 


tree , 


then 


B IS upper 


triangular. 


If G^ 


is a 



£ 

is a one-tree, then b is upper triangular except tor one 
element below the diagonal. Each component can be viewed as 
having one cycle if we assert that each rooted tree has a selt- 
cycle corresponding to its root node. 

In our example basis, the subgraph G^ has only one com - 
ponen t , (one-tree) shown in Figure 2: G^ is a one-tree, with 

nodes 2, 3 and 11 composing its cycle. 

As in the pure network case, this near tr iany u lat ion and 
associated sub-graph are naturally represented by a pre de ces sor 
^nctioji p( ), and a predecessor graph (which does not preserve 
the orientation ot arcs in the original network). The predeces- 
sor function can be used iteratively to construct the unique 
bac kpath from any node to the root (or cycle); the backpath 
includes all nodes on the cycle. The immedia te successors ot a 
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A Preorder Near-Triangulation 
(The underlined coefficient is below the diagonal.) 
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Figure 2. A Generalized Transshipment Basis 
(For the problem in Figure 1.) 
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node, if any, are the first nodes encountered on all paths ex- 
cept the backpath to the root, and all the nodes on these paths 
are called successors . 

Note that each basis may have many near tr ianyu lat ions . 
However, all such near trianyulations yield the same predecessor 
function and yraph (where the riyht to left orderiny of succes- 
sors of any node is immaterial). Thus, the predecessor yraph 
does not completely represent a near tr ianyulation v/ithout 
additional information: an orderiny of the rows (nodes). for 

algebraic reasons, we restrict such partial orderiny to preorder 
[6], in which a node i always precedes its successors, it any 
and in which all its successors, if any, precede any node which 
does not precede node i . 
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IMPLtMtilMTATiUN 



4 . 

For didactic reasons, we begin by introducing a complete 
primal generalized network algorithm using a preorder traversal 
method. Controversial alternatives are deferred until this 
paradigm is presented. Hereafter, notation with upper case 
roman letters followed by parentheses indicates a program data 
array. For instance, the predecessor array is referred to as 
P( ) . 

Static arc storage is used for tails T( ), heads H( ), 
costs C( ), multipliers MUL( ), and capacities CP( ). Contigu- 
ous storage by tail, or by head node reduces T( ) , or H ( ) to 

an hierarchical node-length entry point array. (GENNET uses 
contiguous storage by head node, as does GNET.) Lower bounds 
on arc flow are translated out prior to solution, with appropri- 
ate adjustment of the initial right-hand side of (GT), and of 
CP( ). The sign bit of CP( ) is available to indicate arcs 
nonbas ic at their upper bounds ( ref lected with flow -CP( ) ). 

The predecessor function and its array P( ) are defined 
so that the basic arcs in a cycle are oriented uniformly in a 
directed cycle. All basic arcs not on a cycle are oriented so 
that a backpath is created to a cycle. To obtain this orienta- 
tion the direction of some arcs must be reversed, and the sign 
bit of the predecessor array is used to indicate: if P(I) < 0, 

then the orientation of arc (I, -P(I)) is reversea from its 
original orientation.* 

This is the complement of the discipline used in GNET [6J . 
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A depth array D( ) reveals for each node the number of 



nodes on the backpath before encountering a cycle. Nodes on a 
cycle have depth zero. Number of successors , or preorder 
distance are acceptable substitutes tor depth (6J, but are not 
discussed here. 

A preorder traversal array IT( ) is maintained so that all 
preorder successors of a cycle node are encountered before 
another cycle node. It is convenient to make this a circular 
list for each near triangular basis component by setting IT( ) 
of the last preorder node in the component equal to the first 
preorder node in the component. 

The components of G are not inter-connected, or equiva- 

— D 

lently, the sub-matrices B in (1) do not have common rows or 
columns. Consequently, the p components of a basis may be 
represented in a single set of node-length arrays. 

The array X( ) contains the values of basic variables, 
valups of dual variables (or simplex multipliers, or node poten- 
tials) are stored in U( ), and IVAR( ) gives the location of 
basic variables in the arc arrays. The array FAC ( ) contains 

the cycle factors , defined later. Figure 3 shows these arrays 
tor the basis given in Figure 2 . 

Generalized networks do not exhibit totally unimodular 
bases. Consequently, floating point representation is required 
for X( ), U( ) and FAC ( ), and is desirable for arc-length arrays 

C( ) , MUL( ) , and CP( ) . 
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Figure 3. GENNFT Basis Representation Arrays 
(for basis in Figure 2) 
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step Si, Priceout 



to 



The reduced cost for nonbasic arc k , oriented from fj^ 
t|^ is (given the current dual solution u and column N ); 



''k ’ <=k - 



* -k“t. 
k k 



= C(k) - U(fj^) + MUL(k )*U( t|^ ) . 



(If CP(k) < U, arc k is reflected and the sign of r, is 

K 

reversed.) At most, one multiplication, addition and subtrac- 
tion are required. Note that the multiplication is unneccessary 
I— k^ ~ specialization is possible for sets of 

priced arcs with common attributes. If arc k is a logical arc 
(slack, artificial, or surplus variable) then C(k) can be 
logically generated, rather than explicitly stored, and 



r^^ = C(k) ± U(f|^) , 

depending upon the sign of P(tj^) . 

From the example, 

r^^ - C(Z1) - U(S) + MUL( 27)*U( IJ ) 

= 2U.67 - (27.S52) + .8b(-13.y22) 
= -19.133 . 
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This variable will be used as the entering non-basic variable 
for further illustration. 



Step S2, Ratio Test 

For the determination of the arc to leave the basis the 
system of equations 



BZ 



k 



f 



must be solved for the transformed column Z 
used if arc k is reflected.) Due to the near 
of B , this incoming column transformation can 
with the ratio test in a single integrated proce 
Suppose that N has two nonzero coefficie 
an arc oriented from fj^ to tj^ , and that the 
reflected. An apparent complication arises if 

3 

in separate components of B in (1), say B a 
tively. In this case two disjoint subsystems mu 
the results added to determine the nonzero eleme 
The subsystems are: 



/ K.k 

( -N IS 
triangularity 
be combined 
ss . 

nts representing 

arc is not 

f, and t, are 
k k 

nd B*" , respec- 
st be solved and 
nts in Z 



B^g 




and 






~-k®t, 



18 



k, k. 

(e is the t. th unit vector; y and y are disjoint 

^ k 

components of Z . ) 

This complication is inconsequential. In order to see 
this/ consider solving one of the systems: 



B*'y 



"51k®t. 



The only nonzero elements of y will be those that correspond 

to the nodes in the backpath from node tj^ . As we shall see, 

this follows from the manner in which the coefficient -m, in 

row t propagates during substitution solution. 

K 

Suppose that is a rooted tree. The backpath from node 

t|^ can be denoted by iteration of the predecessor function: t^^, 

P(t. ), p (p ( t, )),... , root; this sequence is shown below as d, 
k K 

d-l,...,U, analogous to the backpath length remaining to the root. 
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The values of a and b for each basic arc depend upon the 
original orientation of the arc, given by the sign of P( ). 
For original orientation, P( ) > 0 implies that a = +1 and b 
is the multiplier value, P( ) < 0 reverses these definitions. 



The triangular system corresponding to this backpath is: 



aQ bj^ 




0 


^1 




0 


• 

• 

• 


c 

ft 

II 


• 

• 

• 


^d-2 *^d-l 




u 






0 






-51k 


• 

• 




u 


• 







By back substitution, its solution is: 



- A 

^d a , ' 



Qf ^6+l*^6+l £ r ^ 1 n 

^6 = for 6 = d-l,...,0 . 

^6 



Now suppose that G*" is a one-tree. Let node tj^ be on the 
cycle. The backpath is tj^ , p(tj^) , p(p(t|^)), ..., c, with 

p(c) = t|^ and length p ; this sequence is shown below as 
-1 , -2 , . . . , -p , analogous to the backpath length beyond the cycle 
start . 
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The corresponding near triangular system is: 







0 


®-p+l 




u 


• 

• 

• 


Q = 


• 

• 

• 


a o b ~ 




0 


”3 ”2 






a ~ b , 




0 


“2 -1 






b a 1 




-m, 


-p -1 




— k 


• 

• 




0 


• 







By modified back substitution, its solution is 



~-k 1 






"^ 6+1 ^ 6+1 



for 6 "2,***, p. 



where 



-p -b. 



f = 1 - n 



6=-l “6 
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The cycle factor (or loop factor [13]), f , is common to all 
nodes in the cycle and can be computed when the cycle is 
created and stored in FAC( ) for all nodes on the cycle so 
that it is immediately available at this step. 

Suppose that the entering arc is (f|^,t|^) with coefficient 
entries (aj^,b|^) and that the fj^ and tj^ backpaths converge. 
Using nomenclature for the backpath sequences introduced above, 
the corresponding system is: 







u 


*1 

• 






• 

a„ , b„ b b 






Jl-1 £ w r 




• 


a ^ 




• 


1 

• 






• 

• 

^w 


^k 

0 ^ 




• 

• 






• 




0 








• 

• 




0 


a 




• 


r 

• 

• 




• 

• 


• 

a j T b , 




0 


d-1 d 






^d 




*^k 


• 

• 




u 


• 




1 
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by back substitution, its solution is: 







(d begins tj^ backpath sequence), 



-b 



(5 + 1^6 + l 



tor 



6 = d , d-1 , . 



^d 

^6 



— , (d begins t. backpath sequence) , 
^d 



^6 + l '^6 + l 



tor 6 = d , d-1 , . . . , w , 



y«-i 



Mr* 

^£-1 



^6 



^(S + l'^6 + 1 



tor 6 = i- A , i- J , . . . , 0 . 



Suppose that entering arc (tj^,t|^) with coetticients 
(aj^,bj^) enters and that the backpaths converge on a cycle, 
corresponding system is: 



The 
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0 



a 

”P 



-p+1 



a 



•p+1 




b 

w 



b 

”P 




-2 




b 



r 



a 

w 





] 

i 



"d-1 



b 

a 



d 

d 



a 

r 



rfl 




0 

0 



Vi ^d 



0 

b. 

0 
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I 



by modified back substitution, its solution is: 



q . = — , (d beyins t. backpath sequence) 

a ^d “ 



^6 = 



^ 6+1 ^ 6+1 



for 6 = d, d-1 , . . . , r , 



^-s-1 = 



-b q , 
r^r 1 

X — . 

a . 1: ' 

-s-1 



*^6 + 1 ^6 + 1 



for 6 = -s-2, ... ,-s , 



q. = — , (d beyins tj^ backpath sequence) , 
^d 



^^6 = 



^6+i ^6+1 



for 6 =d,...,w , 



-b q 

U7 ^ 



‘St-l " i— ^ ^ 



w^w 1 

f ' 



-t-1 



*^6 + 1 ^6 + 1 ^ 

q- = for 6 = 



p, S' , . . . , ^ , 



tor 6 = -i,...,-p . 
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b'or the cycle in Figure 2, 




X 



-( 1) -(-1.17) 

-.79 ^ 1 



4B1U1 



By now it should be apparent that one composite back sub- 
stitution scheme will suftice tor all cases. The cycle factor 
is applied once when, and if, a cycle is encountered on a 
backpath . 

If the backpaths of fj^ and t^^ converge, let joi n be the 
first node on the tj^ backpath that is also on the backpath. 

If the backpaths converge on a cycle and it the leaving arc pre- 
ceeds the join on the fj^ backpath, define join as the first 
cycle node encountered on the fj^ backpath. If the backpaths do 
not converge, join = ^ . 

Several schemes are available for identifying the join 
efficiently [bj. The depth (or number of successors, or pre- 
order distance) of nodes on the backpaths can be used to avoid 
iterating either backpath past the join. Depth, the number of 
nodes on the backpath until a cycle is encountered, can be used 
to indicate which backpath node is deeper and should be iterated. 
When both backpath nodes have matching depths, zero depths indi- 
cate that each backpath is on a root cycle. By remembering the 
first root cycle node on each backpath, further iteration will 
either reveal the root cycles to be distinct with join = $ , or 
coincident with join defined as above. When both backpath nodes 
have matching depths greater than zero, the nodes are compared for 
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equality. A match indicates the join, and a mismatch indicates 
that both backpaths should be iterated for another comparison. 

If a join is encountered, the backpaths have merged and 
either one can be used to complete the ratio test (GENNET 
continues the tj^ backpath). When a cycle is encountered the 
q values are computed on the first pass around the cycle. On 
the second pass around the cycle the q values are computed ana 
the ratio test is completed. 

As the backpaths are iterated, the column transformation 

1 ^ 

is applied and the resulting terms of transformed column, Z , 
are use(;3 in the ratio tes t , seeking the minimum ratio : 

f 

CP(k) the capacity of the incoming arc, 

> 0, node I, 



If a zero ratio is encountered during this process, the ratio 
test may be preemptively terminated. 

Step S3, Pivot 

IP CP(k) is selected as the minimum ratio, then the enter- 
ing variable remains nonbasic and is reflected to its opposite 
bound. Only the flows X( ) need be updated. To do this, the 
backpaths are iterated again ana tor each node l encountered, 

X( ) is reduced by z^ x CP(k) . 



MIN 



X( ^ ) 



for 



CP( IVAK(£ ) )-X()l ) 



for 



- z . 
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If a basis exchange is required, an efficient update of 
the basis representation must preserve the rooted cycle orien- 
tation, updating some entries in P( ) and U( ). Also, some 
flows X( ), and some dual variables U( ), must be changed. 

FAC ( ) must be established for nodes on newly created cycles. 

Some bookkeeping in IVAR( ) and CP( ) may also be required. 

The apparent intricacy of our task is deceiving. Careful 
analysis yields an elegant solution. However, the supporting 
arguments require close attention. 

To simplify the explanation, reorient the incoming arc 
(fk,tk) to (i,j) or (j,i) (if necessary) so that the minimum 
ratio is on the j backpath . Let the entering arc (i,j) have 
the outgoing arc (c,d) on its j-backpath. Also, reorient the 
outgoing arc (if necessary) so that the first node encountered 
on th^ j backpath is c . Figure 2 shows a case for which 
both reorientations are necessary. 

In our example, arc 20, oriented from node 2 to node 11, 
leaves the basis and arc 27, oriented from node 5 to node 13 
enters. We call the backpath segment from j to arc (c,d) the 
j - stem . In Figure 2, i, j, c, and d are shown. The j-stem is 

composed of nodes 5, 2, 3, and 11. The skeletal update: 

a. Reverses the orientation of arcs within the 
j-stem; and 

b. orients the entering arc so that it precedes this 
redirected path with the same orientation. 

If the i and j backpaths merge, the node where they 
merge is called the join. The join is node 3 in Figure 2. 
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If the leaving arc lies beyond the join on the backpaths 
a new cycle is created in the basis exchange. In this case 
the portion of the i backpath from node i to the node with 
the join as its predecessor is called the i-stem. When node i 
is on the j backpath the i-stem is null. In Figure 2, the 
i-stem is node 13. 

Figure 4 displays the pivot logic to be applied. The 
algorithm visits each node affected by the basis update 
exact ly once . It proceeds up the j-stem one node at a time 
visiting the preorder successors of each stem node via IT( ). 

If a join is encountered it switches to proceed up the i-stem 
one node at a time and then returns to the j-stem. At each 
j-stem node, the successors of the next lower stem node have 
already been visited. The unvisited successors of the current 
stem node can be divided into two groups: the left successors 
afe the nodes visited in preorder by iterating IT( ) from the 
current stem node until the next lower stem node is encountered, 
and the right successors of the stem node are the remaining 
unvisited nodes reached by further iteration of IT( ). In Figure 
2, nodes B, 14, 12, and 7 are right successors of node 2 and 
nodes ID, 13, 1, 6, 15, 9, and 4 are left successors of node 3. 

As we climb the i-stem the traversal IT( ) is modified so that 
the last of the left successors (if any) points to the first of 
the right successors and the last of the right successors (if 
any) points to the previous stem node. As we climb the j-stem, 
the traversal is modified so that the last of the left succes- 
sors (if any) points to the first of the right successors and 
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Figure 4. Pivot Traversal Scheme 
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the last of the right successors (if any) points to the next 
node up the stem (because the update reverses predecessor 
orientation for the j-stem). 

The immediate switch to the i-stem upon encountering the 
join during the j-stem iteration is motivated by a subtle 
complication: while visiting the left and right successors 

of the j-stem nodes, the nodes on the i-stem and their succes- 
ors must be skipped if encountered. Because the i-stem nodes 
are successors of the join, visiting the i-stem as soon as the 
join is encountered (if one exists) on the j-stem leaves us 
with the preorder successor of the last i-stem node visited. 

This valuable artifact enables the subsequent j-stem iteration 
to immediately skip all i-stem nodes and their successors should 
they be encountered. This is the key step preserving an effi- 
cient one-pass basis update. In Figure 2, i-stem node 13 and 
its successor node 1 are successors of the join, node 3. The 
preorder successor of the last i-stem node (called the preorde r 
link in Figure 4) is node b. 

A stem node may have a right successor which is on the 
root cycle (with depth zero). The preorder traversal array is 
organized so that all successors of a cycle node are encountered 
before the next cycle node. This implies that a cycle being 
broken by a leaving arc will always be encountered as a right 
successor of a stem node. In Figure 2, the broken cycle is 
encountered as a right successor of node 2; if arc (2,11) were 
not the leaving arc, then node 11 would be a preorder successor 
of node 2 . 
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The basic arc flows, X( ), are changed as each arc is 
visited on the j-stem, and (if a new cycle is created) on the 
i-stem. If no cycle is created, X( ) is changed only it the 
minimum ratio is nonzero, and then only on the arcs visited on 
the backpaths. 

During the update, the dual variables must be changed so 

that 




for every basic arc k oriented from node tj^ to node f^^ . 
With incoming arc (i,j) (reoriented as in Figure 2 so that the 
outgoing arc is on the j-stem), this relationship is retained 
for all nodes except for those which the update changes to be 
successors of i: (i.e. the nodes on the j-stem and their 

successors ) . 

If a cycle is not created, this update proceeds for eacii 
j-stem node and its successors as these nodes are visited in 
preorder. For node s , and associated basic arc £ oriented 
from s to p(s) , 

U(s) = C(£) + MUL(£) * U(P(s)) , 

while the reverse orientation -P(s) to s yields 
U(s) = (C(£) - U(-P(s) ) )/(-MUL(£) ) . 
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As in pure networks, the preorder traversal assures that a value 
of the dual variable of a predecessor node is always determined 
prior to its use by any immediate successor nodes. 

When a cycle is created, the dual variable must be deter- 
mined for one of the cycle nodes. Then, the dual variables of 
the remaining cycle nodes and their descendents can be found 
one at a time in preorder traversal. 

This key dual variable is computed immediately after the 
ratio test predicts creation of a new cycle (the leaving arc 
lies beyond the join on the ratio backpaths). Consider the 
current context for a new cycle: 



/'^N Join 




from which the new cycle will be formed (with modified pre- 
decessor function p( )): 




p(p(j)) P(j) 
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The corresponding near-triangular system is: 



(u 



-p 



u 



-p + 1 



a b , , 

-p -p+1 



*-p + l 



, u_ ) 



“>-2 





/ C 



-p + 1 



'-1 



) 



where the indexing of c is understood to yield basis arc 
costs (accomplished by indexing with IVAR{ ) ). 

The determination of one term (say, ) of the solution 
of this system is induced by modified forward substitution to 
be ( e ,g . , [ 11 ] ) ; 




where (the cycle factor) f is defined (again) as 




and u_|^^ is the corresponding term of the solution of the 
upper triangular system (omitting the cyclic coefficient 
and using forward substitution): 



strictly 



b 

-P 



f 
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c 



u 



-p a 



u 



6 




for 6 = -p+1 



1 



Note that computation of requires that we traverse the 

new cycle in a direction oppos i te to its predecessor orientation 
However, before the update creating the new cycle, the j-stem 
exhibits proper orientation for at least part of our work. Thus 
we can complete the first portion of the forward substitution 
for and accumulate the associated partial product coiaponent 

of f while iterating the j-stem before the update. 

The remainder of the new cycle is accessed by iterating 
the i-stem. As we proceed up the i-stem the remaining product 
terms of f are accumulated, and the reverse i-stem path is 
stored (e.g., using U( ) locations, which contain obsolete dual 
values to be replaced during the imminent update). Reaching 
the join, this stored reverse i-stem path is then accessed to 
complete computation of uj^^ . 

The newly created cycle in Figure 5 is composed of j-stem 
nodes 5, 2, and 3, and i-stem node 13. The dual system becomes 




4b . 6U 




lb .47 



u 



2 



74u 



5 



26.76 



88u 



u 



20.67 



13 



5 
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which yields 



= -19.0142, 

f = 0,3488 (the new cycle factor), and 
~ -54.5132 (the new dual for node 13). 

Thus, when a new cycle is to be formed, the new cycle nodes 
must be visited once (after the ratio test and prior to the 
pivotal update). 

The one-pass preorder traversal update can now proceed as 
presented in Figure 4. The basis representation arrays are all 
modified on-the-fly during this traversal. The update of nodes 
on a newly created cycle (if any) includes establishing the new 
cycle factor FAC ( ). Changes for P( ), D( ), IT( ), IVAR( ), 

X( ), and U( ) proceed analogously to the pure network case 
(e.g., GNET [6]) with simple modification of generators for X( ) 
and U( ) to accommodate generalized network coefficients for 
basic arcs. 

Figure 5 shows the new basis (derived from the example in 
Figure 2) before restoring near-triangulation with the update. 

A new cycle is to be created and the new cycle factor and a dual 
solution (for node 13 on the i-stem) are found at this stage. 

Figure 6 displays the new basis restored to near-triangu lar 
form. At this point, all basis representation arrays are updated 
with the values shown in Figure 7 (data in italics has been 
changed by the update ) . 
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(The entering column is ltallclze.cL) 









Join 




Figure 5. New Basis Before Restoring Near-Triangulation 
(entering column 27, arc (5,13)) 
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A Preorder Near-Triangulation 




Figure 6. New Basis Near-Triangulated 
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GhlNNET is designed to exploit intrinsic pure network struc- 
ture commonly found embedded in generalized network problems. 

Note that when this algorithm is applied to pure network prob- 
lems, automatically adapts to a minor variant of GNET. ut 
course, GENNET uses floating point arithmetic operations which 
are intrinsically slower on most computers than the pure additive 
integer arithmetic of GNET (also, floating point arithmetic 
requires some extra editing for mantissa truncation errors). 

To mitigate this disadvantage, GENNET can test logically for 
pure network arcs (with unit multipliers) and avoid unnecessary 
floating point multiplication and division operations. 

GENNET also employs an automatic dual basis aggregation 
refinement ([6],p.26 ff). Explicit values for b( ), U( ) and 
IT( ) are maintained only for nodes with successors. An’ array, 

A( ), records for each node the current number of its aggregated 
successors . When an aggregated node is encountered in the 
priceout, its dual is generated from that of the immediate 
predecessor of the node. When a backpath of an entering arc 
begins with an aggregated node, it is disaggregated, and when the 
leaving arc isolates nodes with no successors, they are 
aggregated . 

Figure B shows the arrays affected by the dual aggregation 
scheme for the basis in Figures 2 and J. TT( ) indicates aggre- 
gated nodes with U entries, and these nodes have broken outlines 
in the one-tree depiction. The priceout of arc (5,li) requires 
that U(5) be generated using the predecessor dual U(2). Node b 
subsequently starts the j-backpath and is disaggregated. The 
update leaves node 11 with no successors, and thus aggregated. 
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Figure 7. GENNET Basis Representation Arrays 
(for Basis in Figure 5) 
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Figure 8 . Aggregated Basis Representation 
(for Figures 2 and 3) 
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COMPUTATIONAL LXPLKILNCL 



b . 



Significant design alternatives for GLNNLT have each been 
evaluated by extensive experimentation at large scale. Illus- 
trative computational experience is abstracted in this section 
tor some of the prototype systems tested. Departing somewhat 
from the style of the paper documenting such work for GNBT lb], 
relative performance is reported even tor some competitive design 
features subsequently rejected tor adoption (some readers of (bj 
have concluded, quite incorrectly, that only those features 
reported tor GNLT were tested). This should help other re- 
searchers avoid our mistakes, and may even change some widely 
held misconceptions and correct a few translation errors in 
textbooks . 

Among the key issues to be resolved are: 

a) Static Storage of Arcs . The arc lists can be stored 
in arbitrary sequence, or, to save space, arcs can 
be stored contiguously by tail node, or by head 
node, thereby replacing an arc-length index array by 
a (presumably much shorter) hierarchical node-length 
entry point array. 

b) Preorder Manipu lation of Basis . The triple label (aug- 
mented predecessor index) method [ZU,22J will be 
presented and compared with the preorder traversal 
method . 

c) basis Aggregation . An aggregated basis representation 
will compete with an explicit representation. 
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d) Pric ing Sch emes . Candidate list schemes and explicit 
arc pricing mechanisms widely used in general linear 
programming systems will vie with dynamic candidate 
queue disciplines. 

e) Pure Networ k Speciali zat ion . Generalized network 
algorithms would ideally adapt to pure networks with 
efticienqy comparable to pure network codes. 

t) Star ting , Tun i n g and Tailoring. Which algorithm 

parameters and settings lead to high elticiency for 
interesting classes of problems? Are heuristics tor 
advanced starting solutions worthwhile? 

g) Gen erali zati ons. Advanced features and generalizations 
will be suggested. 

Computational tests have been made with many problems, 
including a benchmark suite of pure network problems generated 
with NETGSN [28J, and generalized network problems generated 
by NETGENG [17,lbJ. figure y gives some problem characteristics. 

Static arc storage has been implemented in three ways: 

• Contiguous by head node with hierarchical node- 

length entry point array, 

• Contiguous by tail node, with hierarchical node- 

length entry point array, and 

• Explicit arcs in arbitrary order. 

In competition with our basis manipulation using preorder 
traversal, a triple label representation originally suggested 
by E. Johnson [2fc>j has been implemented tor pure networks and 
called tiie augmented predecessor indexing method by Glover, 

Karney , and Klingman 12UJ. Glover, Klingman and Stutz 122J 
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t'iyure y. bpme Benchmark Problems 

(Problem NG27 is extensively studied in [bj.) 
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^report that the method has been extended to yeneraiized networks, 
but reveal no details. 

We have implemented our own etticient version of the triple 
label scheme. 

The triple label representation uses preoecessor, successor, 
and brother pointers tor each node. figure lU shows these arrays 
tor the basis in figure 2. 

To briefly illustrate the triple label scheme tor gener- 
alized networks, let our situation be the same as tor the pre- 
order example (the j-backpath of the incoming arc is arranged 
to include the outgoing arc and to encounter node c first on 
that backpath). Let 




r 



and the backpath from j to c 



v 



d' d-1 



c . 



The skeletal upoate scheme is: 
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Figure lu. Triple Label Basis Representation 
(for Figure 2) 
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1. Set 6 = d 

2. It ^ t o Ste p ^ 

Otherwise, it possible, tind a node v* 
such that P(v*) = v^_j^ 
and B(v*) = V., 

and set B(v*) + B(v^), Oo 1 

J. Set f- B(v^), 

4. set P(v^) ^ 

B(v^) - 

b. Set 6 + 6-1, 

Then, it 9t c, ^ 2, 

Otherwise, Stop . 

(Step 2 exhibits the key extension tor generalized networks of 
the pure network triple label scheme.) 

These tive steps reverse the orientation ot all arcs on 
the j-stem and orient ^ ^ so that it begins this redirected 
path. All other triple label operations are obvious alterations 
ot the preorder traversal procedure. 

A static candidate list pricing strategy (e.g., [i7,iB,bl])) 

an explic it arc pricing method reminiscent ot general linear 
programming systems, and a dynamic candi date queue lb] have been 
tested . 

The (Lj^,L^) candidate list procedure is a simple strategy. 
Nodes with leaving arcs (or entering arcs for contiguous head 
node arc lists) are sequentially priced, placing the most nega- 
tive candidate (it any) from each node on the candidate list 
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until candidates have been located (arbitrarily organized arc 
lists require explicit arc pricing). t;ntering arcs are chosen 
trom the candidate list by most negative reduced cost until 
iterations have been carried out, or until all list entries have 
non-negative reducea costs. Arcs with non-negative reduced costs 
are dropped from the list and replaced by continuing to scan the 
nodes until candidates have been tound. It an exhaustive pass 
through the nodes results in less than candidates, then an 
optional closing gambit sets equal to the actual candidates 
found, and reduces by halt unless equals one. Glover, 
Hultz, Klingman and Stutz [17,18J report (L^,L^) of (b,lU) to be 
best in their work. 

The candidate queue is a dynamic list of interesting arcs 
and nodes, scanned in a cyclic manner. The entering arc is 
selected from the queue by pricing NNK entries; if an interesting 
node is encountered it is replaced by its best-priced entering 
arc (or leaving arc tor contiguous tail node arc lists). Arcs 
pricing favorably are retained in the queue. When the end of 
queue is encountered, the queue is refreshed by pricing IPG nodes 
in a cyclic general arc scan. During an opening gambit of IMNS 
pivots, the nodes incident to the entering basic arc are added to 
the queue. There is no closing gambit, since the queue automa- 
tically shrinks and finally collapses at optimality. Bradley, 
Brown, and Graves [6] suggest NNh) = 32, WNb = 3m/4 and 
IPG = m/lU + I for pure network problems with m nodes. 

A rule to break ties in the ratio test which guarantees 
finite convergence for pure transshipment problems has been 
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developed by Cunningham [12], Bradley, Brown and Graves 16] 
show that the conditions necessary for finite convergence are 
naturally satisfied by GNBT on over 90% of its degenerate pivots. 
Elam, Glover and Klingman [lb] have observed that the results of 
Cunningham can only be extended to the generalized network case 
when the multipliers are positive. We have not used Cunningham's 
modification . 

bland [4] presents a class ot restrictions of pricing and 
ratio tests for general linear programs which relies exclusively 
on primal simplex representation and guarantees finite conver- 
gence. These rules are easily modified to produce an efficient 
finite simplex algorithm. The modification interferes with ef- 
fective pricing strategies only during degenerate pivot sequences 
and the restrictions increase in severity only with the number of 
pivots in that sequence. However, during a degenerate pivot se- 
quence restriction records must be accumulated (e.g., a list with 
each incoming variable in one ot our schemes). This record is 
naturally accommodated by the dynamic candidate queue, but not by 
a static candidate list or explicit pricing. No purpose is 
served by reporting such unbalanced competition. 

A starting strategy has also been tested in conjunction with 
pricing alternatives. A straightforward starting method ex- 
amines each node with supply (or demand tor contiguous arc stor- 
age by head node) and assigns as much flow as possible to its 
least cost leaving (entering) arc. The procedure stops when an 
exhaustive pass of the nodes makes no additional flow assignments 
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The starting solution achieved is not necessarily teasible. 

(E.y., "exhaustive pass sequential source minimum start" [laj.) 

Artificial arcs are driven from the basis in all experiments 
using a Big-M method (e.g., [b]). This choice is principally 

motivated by the comparability of competitive tests between pure 
and generalized network codes on pure and generalized network 
problems. (A two-phase method is employed in production use.) 

Choosing the best Big-M value is a bit tricky. The smallest 
Biy-M value which yields a feasible optimal solution (it one 
exists) is best in our experience. Bmall Big-M values may tail 
to produce feasibility, and large values inflict numerical 
difficulties. In practice, a default value is used and an 
automatic restart recovery is applied it an infeasible solution 
persists. It a restart with a higher Big-M value tails to 
reduce the total infeasibility, a terminally infeasible solution 
is declared. Figures 11 and 12 indicate the multiple of maximum 
absolute arc cost used for Big-M in each problem. 

Computational tests have been performed on various computer 
systems. The times reported here are accurate to the precision 
displayed tor IBM 37u/lby-B using the FORTRAN H compiler with 
UPT(2). Solution times exclude input/output overhead. OBNNBT 
uses double precision (IBM RBAL*y) arithmetic tor floating point 
operations and storage. 

Solution times are given in Figure 11 for the pure network 
test problems. Performance is given tor three pure network codes 
(two versions of ONBT and suPBRk) as well as tor several repre- 
sentative generalized network prototype systems. we can thus 
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compare the best generalized network scheme with the best pure 
network code exhibiting equivalent features (GNbT, or its 
aggregated successor variant [6J). 

The times tor SUFtlRK also provide the only available objec- 
tive means tor comparison of our implementations to that ot 
Glover, Hultz, Klingman and stutz [14]; they report that their 
generalized network code, NbTG, is about as fast as sUFtKK (a 
fast out-ot-kilter code for pure networks (Ij) when both are used 
to solve pure networks. Using their version ot yUPUKK on our 
computer we have shown that our implementation of NUTG (called 
TLA in our nomenclature) is at least as efficient as their claims 
for NETG. 
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Figure 11. Pure Network Test Problems 
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t ifjure 1.:. irar^s lipment (G'l ) Test Problems 



However, in these tests TLA is substantially outperformed by 
the alternate systems. For the pure network problems, best 
performance is achieved by (Figure 11); 

Contiguous arcs by head node. 

Candidate queue pricing. 

Preorder traversal, and 
Aggregated successors. 

TLA does not incorporate any of these features, and is generally 
less than half as fast as competitors. 

Although GENNET (HUPX) should in theory rival GNET (HgPX) 
with pure network problems, the overhead of testing in GENNET 
for more general basis structure and the additional computa- 
tional burden of floating point arithmetic exact a performance 
penalty of about 15 percent. 

Figure 12 shows solution times for the generalized trans- 
shipment network problems. Note that the starting strategy helps 
candidate list performance and hinders the candidate queue. 
Arranging arcs contiguously by head node dominates both tail 
node and explicit arc list designs. The candidate queue pro- 
vides good performance accompanied with contiguous arcs by 
head node. Preorder traversal continues to provide better 
performance than triple label representation in all design 
contexts. Aggregated successors offer a pronounced advantage. 

GENNET (HUPX) provides best overall performance. It 
offers a decided advantage on problems with many more sinks 
than sources, a situation common in real life. 
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Figure 13 displays performance of GENNFT(GN, HPQX) applied 
to a set of (GN) problems extracted from a collection of real- 
world LP/MIP models [lUJ. Despite the slight additional floating 
point arithmetic required to solve (Gb), GENNET solves these 



problems much faster 


than would 


be pred 


icted by experience with 


randomly generated GT 


problems . 


Tuning 


of the pricing mechanis 


greatly enhances this 


difference 


• 






Problem 


Nodes 


Arcs 


Seconds 


Pivots 


AIRLP 


17U 


3 ,040 


2.62 


420 


COAL 


17U 


3,923 


1 .80 


471 


STEEL 


422 


1,279 


.39 


499 


FOAM 


y5i 


4,953 


3.74 


1 , 258 


ODSAS (GN ) 


1,431 


4,615 


3.22 


1,427 


ALUMINUM(GN ) 


2,17B 


7,216 


3.57 


2,794 


REFINE (GN ) 


3,110 


6,617 


4.72 


3,322 


FOOD(GN ) 


3,716 


13,907 


12.11 


7,004 



(Big-M = lU X largest cost coefficient) 



Figure 13. (GN) Test Problems 
(GN rows extracted from real-world LP/MIP models [lU] 
with null columns deleted and slack arcs added.) 
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Close scrutiny of solution trajectories lends some insight 
into GENNET's good performance. GENNET has a one-pass itera- 
tion unless a cycle is formed; a cycle is formed on only b-to-24 
percent of all iterations for these problem sets — 5-to-lU percent 
for most problems. Also, the explicit ( non-aggregated ) subset of 
the nodes is remarkably small, seldom numbering much more than 
the number of source nodes. Finally, the length of backpaths is 
quite short, averaging about the number of echelons (path length 
from sources to sinks) in the model, or just more than 2 in these 
problems . 
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6 . 



CONCLUSIUN 



The generalized network system Gh'lMNt;!' is small, tast and 
easy to modity. Adaptations have already included using GtNNhT 
in a system to solve generalized networks wiuh complicating side 
constraints an/or complicating variables (McBride [3Uj). GBNNt'i' 
has also been incorporated in a powerful microcomputer-based 
network optimization system by brown, Dutt and t'inley [7,16J 
using an APPLE-II host and PASCAL implementation language. 
Modifications for mixed integer generalized networks have also 
been tested (though not with care sufficient to warrant 
publication at this time). GLNNET has proven to be a worthy 
successor of GNET [6]. 

Preorder traversal is appealing for its mathematical and 
implementation elegance, and has proven to be efficient and 
flexible for generalized networks (as it was tor pure networks). 
(Adolphson and Heum [1] have also suspected this and have in- 
dependently pursued this avenue.) 

Experience shows that the Gennet design performs much more 
efficiently on real models than on randomly generated test 
problems of nominally equivalent size; this design is also tech- 
nically and philosophically compatible with the various systems 
we have devised tor solving other more general classes of 
optimization models. 
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The FUKTRAN programs GENNET — (GT) and (GN) versions — 
(Copyright iyb4) are licensed to researchers tor a nominal 
charge on an exclusive use basis. For further information 
the authors via P.U. Box 1832, Alexandria, Virginia, 22313, 



write 

USA. 
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